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Abstract 

We find an explicit solution of the Schrodinger equation for a Chern-Simons theory 
coupled to charged particles on a Riemann surface, when the coefficient of the Chern- 
Simons term is a rational number (rather than an integer) and where the total charge is 
zero. We find that the wave functions carry a projective representation of the group of 
large gauge transformations. We also examine the behavior of the wave function under 
braiding operations which interchange particle positions. We find that the representa- 
tion of both the braiding operations and large gauge transformations involve unitary 
matrices which mix the components of the wave function. The set of wave functions are 
expressed in terms of appropriate Jacobi theta functions. The matrices form a finite 
dimensional representation of a particular (well known to mathematicians) version of 
the braid group on the Riemann surface. We find a constraint which relates the charges 
of the particles, q, the coefficient of the Chern-Simons term, k and the genus of the 
manifold, g: q 2 (g — l)/k must be an integer. We discuss a duality between large gauge 
transformations and braiding operations. 
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It is by now well established that particles confined to a two dimensional space can have 
fractional statistics. Such particles are called anyons. Interest in them is partially motivated 
by their physical effects such as their conjectured role in the fractionally quantized Hall 
effect or high temperature superconductivity, and partially by the fact that the description 
of anyons uses interesting mathematical structures. Anyons are a generalization of ordinary 
bosons or fermions where the wave functions of many identical particles, instead of being 
symmetric or antisymmetric, carry a representation of the braid group on the two dimensional 
space. This is often described mathematically by coupling the currents of the particles to 
the gauge field of a Chern-Simons theory, 

S = -^jAdA + J A,fd 3 x (1) 

The solutions of this field theory take the form of a set of wave functions which at a given 
time depend on the particle positions. It is well known that these wave functions carry 
a non-trivial representation of the braid group, where the representation is specified by 
the parameter k as well as the number of particles and the topology of the 2 dimensional 
space. The braid group is an infinite discrete non-Abelian group and has many potentially 
interesting representations (see, for example H). 

Abelian Chern-Simons theory coupled to classical point particles was solved on the plane 
by Dunne, Jackiw and Trugenberger M. In that case there are no degrees of freedom, the 
Hilbert space is one-dimensional and the only quantum state is given by a single unimodular 
complex number. The phase of the complex number depends on the parameters of the 
model, i.e. the coupling constant k and the histories of time evolution of the positions of the 
particles in such a way that if we consider a periodic trajectory, Zi(t), i = 1, . . . , N, t G [0, 1] 
and Zi(l) = Zi(0) the phase of the wave function changes by the well-known factor 




which counts the changes of relative angles of positions of the particles. This can be inter- 
preted as the wave function carrying a one-dimensional unitary representation of the Nth 
order braid group on the plane. The statistics parameter is the change in phase of the wave 
function under interchange of two particle positions. Here it is given by f . 

In another interesting paper, Bos and Nair || have solved the Schrodinger equation for 
Abelian Chern-Simons theory coupled to particles when the space is a Riemann surface of 
genus g and when k, the coefficient of the Chern-Simons term, is an integer. In this case 
the representation of the braid group which is obtained is more interesting. It is generally 
multi-dimensional and represents the braid group on the Riemann surface which has more 
(non-Abelian) structure relations than the braid group on the plane. In this Letter we 
shall show how this solution generalizes to Chern-Simons theory with k a rational number 
and examine the resulting properties of the braid group representations which we obtain, 
realizing explicitly the representations found in [|]]. 

In Abelian Chern-Simons theory, unlike the non-Abelian case, there is no reason why k 
should be quantized but could in principle be any real number. Polychronakos [fj has shown 
that, even on a torus, k need not be quantized. Furthermore, Witten has argued that 
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the dimension of the Hilbert space of pure Chern-Simons theory on a Riemann surface with 
genus g is k 9 when k is an integer. Polykronakos || showed that when k is rational, the 
dimension of the Hilbert space of pure Chern-Simons theory is (k^) 9 - It is also known that, 
when k is an integer, the representations of the braid group on a Riemann surface obtained 
by coupling to point particles depends on k and g [|J. 

We start from Abelian Chern-Simons theory coupled to the current due to a gas of charge 
particles (|IJ). We shall first demonstrate how to extract the physical degrees of freedom 
from the non-trivial cohomology classes of the gauge field and the complex structure of 
the Riemann surface. Later, invariance under large gauge transformations and modular 
transformations will be used to show that the Hilbert space is finite dimensional. 

Before quantizing, we rewrite (|l|) in term of topological quantities on the 2-dimensional 
Riemann surface, A4, of genus g. We decompose A into its exact, coexact and harmonic 
parts. More precisely, the Hodge decomposition of A, on M, is given by (d and * act on M. 
in this paper) 

1 1 ?7T? 9 

A = d(— *d*A) + *d(—*dA) + — £(W - W) (2) 

where I/O' is the inverse of the laplacian (□) acting on 0-forms where the prime means that 
the zero modes are removed. The zero modes of d and *d (or □) acting on a one-form are 
spanned by the set of Abelian differentials, u l , on J\4. 

We can represent the homology of M. in terms of generators a\ and its conjugate gener- 
ators bi, I — 1, . . . , g. The intersection numbers of these generators are given by 

v(a h a m ) = v(b\ b m ) = 0, u{a h b m ) = -u(b m , a t ) = 5? (3) 

where v{C\ 1 C<2) is the signed intersection number (number of right-handed minus number 
of left handed crossings) of the oriented curves C\ and C%. The holomorphic harmonic 
one-forms uj 1 have the standard normalization [ij 

u m = 57 1 , I u m = Q lm 

ai Jb l 

The matrix Q lm is symmetric and its imaginary part is positive definite. This actually defines 
a metric in the space of holomorphic harmonic forms 

if uj 1 A cu m = 2Im(^ m ) = G lm , G lm G mn = 5? (4) 

We will use G[ m and G lm to lower or raise indices when needed and use Einstein summation 
convention over repeated indices. 

Any linear relation, with integer coefficients, of a/ and bi that satisfy @ is another valid 
basis for the homology generators. These transformations form a symmetry of the Chern- 
Simons theory and comprise the modular group, Sp(2g, Z): 




where S = y b A ) ^ 
with SES 1 = E and E = I " " I . The g x g matrices A, B, C, D have integer entries. 
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We can define 

£ = -— — *d*A, F = *dA 
So we then have (including the A term) 

2tt 1 



A = A dt - —d^ + *d(—F) + 2m(^LU l - 7l O l ) (6) 

Similarly we can write the current j = J M gfjr into a one-form j^dx^ = jodt+j, using the metric 
which we assume to be flat in the time direction (goo = 1, goi — go2 = and the remaining 
components forming the metric on M). We can use again the Hodge decomposition of *j 
on M 

*j = -d X + *diP + i(j l D l -j l u l ) (7) 
The continuity equation ^fd 3 x + d*j A dt = can be used to solve for if) 

We shall consider a set of point charges moving on A4, with trajectories Zi(t) and charge 
qi, where Zi(t) ^ Zj(t) for i ^ j. For technical reasons (involving the solution of Gauss' 
law) we shall assume here that the total charge is zero, J2% <li = 0. We believe that this is 
only a technical restriction and can be generalized. We shall discuss this generalization, the 
inclusion of the nonzero chrge sector, in a future publication.^] The current is represented by 

Jo{z, t) = Y, QiK z - z M)i j( z > t) = ^i. z - ^( t ))^( i i( t ) d ^ + z~i{t)dz) (8) 
i i 1 

Integrating (§) with the harmonic forms uo l , we find the topological components of the 
current in ([?]) 

j\t) = Y J q*muj\z l {t)) , f{t) = J2u z i(t)u l ( z i(t)) (9) 

i i 

To solve for x, it is best to use complex notation 

R = if) + ix = R(z, z) 
where we find *dx + dip = d z Rdz + d^Rdz. From (0), (||) and using (|9|) we find 

d z R + jiu\z) = - ^qiii5(z - Zi{t)) 



1 Gauss' Law (see ahead (p^7j)) implies that the (single) additional degree of freedom associated with the 
nonzero charge should commute with those in the zero charge sector. This (and the additive structure of 
the Hamiltonian) would imply that the full wave function is a product of the zero charge wave function (see 
ahead equation ([56])) and an additional factor for the nonzero charge, depending only on the single quantum 
variable £o conjugate to J F, and fixed by Gauss' Law to be exp — ^p<3£o- Thus our results will be largely 
anaffected by this inclusion of the zero charge sector. 
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d- z R + 3l u\z) = ^q i z i 6(z-z i (t)) (10) 

i 

To solve ( [TOD for R, we will need the prime form 

E(z, w) = (h(z)h(v,))-i ' e ( J/2 ) (jf w l n ) 

where = ^r© ^ jy^ ^ ( u l^)|u=o -oj 1 (z). The prime form is antisymetric in the variables 

z and w and behaves like z — w when z « it; (the /i(z) which appear in the denominator are 
for normalization) .0 

The theta functions M are defined by 



where a, [3 € [0, 1], have the following property 



i7r(n ; +o ; )C im (n m +a m )+2m( ni +a t )(z l +/3 1 ) 

i 



(11) 



for integer-valued vectors s m and ti. For a non- integer constant c 

e ^ ^ (z m + cn ml u\n) = e^^u-t^^+WQ ^ a ~ ct \ ( z p) 

The solution of ( |10|) is 

B =s l T& l ^ifel )l " ,,|,) l ( " (12) 

where we have chosen R such that R(zq, z ) = for an arbitrary point zq (We can choose 
Zq = oo for genus zero). The important fact about R is that it is a single- valued function. If 
we move z around any of the homology cycles, R returns to its original value. In fact, this 
is also true for windings of Zq, an important relation since it is only a reference point. So 

x - |l-s?«(fgi>>] + sttftM (13) 

Now we are ready to solve the Chern-Simons theory on Ai. By putting (H) and (0) back 
into (|I]) we find 

S = ~ / (£PF - £PF)d 3 x + ink J ^% - 7'7,)dt + / 4>0'o - 



2 This formalism can also be extended to include the sphere, where there are no harmonic 1-forms at 
all (the space of cohomology generators has dimension zero) by properly defining the prime form. We use 
stereographic projection to map the sphere into the complex plane and use E(z, w) = z — w as the definition 
of the prime form. 
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- / + FP X )d 3 x + 27rz J (ja l - j l 7l )dt (14) 

where P is a projection operator onto the space orthogonal to the zero mode. From this we 
obtain the equal-time commutation relations of the quantum theory 

[£(*), F(w)} = -iPS(z - w) or F{z) = zP^y (15) 

and 

The projection operator in ( [TBI) changes the delta function to S(z — w) — l/Area(.M). The 
functional derivative must also be defined using this projection operator. With this holo- 
morphic polarization || it is convenient to use the following measure in 7 space 

(ttiltfa) = f e -2^ 7 ™ Gmi f ^*( 7 )^ 2 ( 7 )|G|- 1 nd7 m d 7 m 

J m 

where \G\ = det(G mn ). With this measure, we find that 7^ = 7 as it should be. 
Aq is a Lagrange multiplier which enforces the Gauss' law constraint 

F W -|* W =«.^-| AW «0 (17) 

Note that since the total charge is zero, the total magnetic flux should be zero, as is necessary 
on a compact space when the gauge field A is a globally defined 1-form. 

Under a modular transformation, the basis 7', •y 1 will be transformed accordingly. This 
will not change the choice of polarization, since the modular transformations do not mix 7 
and 7. 

From (14), ([15]) and (|T6|), we find that the hamiltonian, in the Aq = gauge, can be 



separated into two commuting parts 

H = - I AA*j =H + H T 
Jm 

where 

*' = I<T«* + * p & d '* (18) 

(note that d 2 x = —^dz A d^) while the additional part that takes care of the topology is 

H T = i(2n] n l - p 1 ^) (19) 

To solve the Schrodinger equation, we will use the fact that the hamiltonian separates, 
thus writing the wave function as 

tf(£,7,f) = *o(£,*)*T(7 5 *) 
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with the Gauss law constraint (IT 



which is solved by 



(tf^-yio)*o&*) = 



= exp[-— / f(z)j CM)d 2 a;]tf c (t) 
ft J jvi 



The first Schrodinger equation is 



dt 

which has the solution 



m \ k dt 



d x 



tf c (t) = exp 



2ni r f 



X {z,t') Jo {z,t')d 2 xdt' 



k Jo Jm 

For a system of point charges, the use of (|i"3f) allows us to write (|20D and (^Tj) as 

2vri 



# (£,t) = exp 



2 Z.J 



where 
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m&t* i ] i (t")dt" 



-mat' I 3 l (t")dt" 



+ 



7T 

2k 



and 



8ij(t) = Imlog 



E(z i (t),z )E(z ,z j (t)) 



+ Im 



2i(0) 



2j (t) 



+ LUi) + 



*j(0) 



0/ 



Zi(t) 

Zi(0) 



(20) 



(21) 



(22) 



(23) 



(24) 



is a multi-valued function defined using the prime form. We will need the phase (|23|) for the 
topological part of the wave function. The function 6^ (i) is the angle function for particle i 
and j. For i = j, we choose a framing Zj(£) = Zj(t) + e/j(£) which lead to the replacement of 
E(zi(t),Zi(t)) by fiif). So the wave function (|22|) , with the angle function (0), accurately 
forms an Abelian representation of the braid group O, 3 . 

Now, the topological part of the hamiltonian is used to find the part of the wave function 
affected by the currents going around the non-trivial loops of M.. The Schrodinger equation 
for (M) is 



dt 



H T ^ T {l,t) = i \ 2irja l - -j 







k 



*t(7,*) 



which has the solution 



*t(7,£) = exp 



27T7' / Ut')dt' - 



2tt /•* 



k Jo 



m&t' i ] i (t")dt" 



*r(7.*) 



(25) 
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Note that with the phase (|23|), the double integral above will turn into /q ji(t')dt' ■ Jq j l (t')dt' , 
a topological expression. 

The remaining equation for ^(7,^ 



dt 



~k 3 



d'y 1 



is easily solved in the form 



i> T { 1 \t) = ^-ljj\t')dt') 



(26) 



(27) 



The wave function (27) is not arbitrary, but must satisfy the invariance of the action (|l|) 
under large gauge transformations, when there is no current around. So let us set = for 
a while and find the condition on ^t- 

In general, the large U(l) gauge transformations are given by the set of single- valued 
gauge functions, with s m and t m integer-valued vectors, 



U s>t (z) = exp (2m(t m ri m (z) - s m r] m (z)) 



where 



v m (z) = % / (ft"% - QTQx) 



fj m (z) 



(COr, 



LU r , 
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If we change the endpoint of integration by z — > z + aiu + b m v m with u, v integer and a, b 
defined in fl3|), we find tf 1 — > 77™ + u m , fj m — > ry m + f m and f/ Sjt — > U Sjt e 2ni< - tmUm ~ smVm ^ = U S)t . 
The transformation of the gauge field (0) under i is given by 



rfi _ 7 m + g m + Q ml ti ? 7™ _> 7™ + s ™ + jym^ 

The classical operator that produces the transformation (|28|) 



(28) 



c s ,t(7,7) = exp 



V U d^ m v ^7" 



must be transformed into the proper quantum operator acting on the wave function ty?- 
By using the commutation ( [TED to replace by —2Txk^ m we find the operators C S)t which 
implement the large gauge transformations || 



C,M = exp -2nk(s m + lrt,)7 m - Kk(s m + n m %)G mn (s n + Q nl t, 



(29) 



The quantum operators C Syt do not commute among themselves for non- integer k. From 
now on we will set k — |r for integer k\ and k 2 . Now, in contrast with their classical 
counterparts, the operators C st satisfy the clock algebra 



L/ S 1 ,t 1 { ^S 2 ,t 2 — e °S 2 ,t2 C 'Si,ti 

Their action on the wave function is 

C s M^T(l m ) = exp \-2vk{s m + n ml U) lm 



(30) 
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-nk(s m + n m %)G mn (s n + n n %)\ ^ T ( 7 m + s m + Q m %) (31) 

On the other hand Ck 2S ,k 2 t commutes with everything and must be represented only by phases 
4> s t- This implies, using fl3T|), 



^r(7 m + h(s m + n ml U)) = exp -i<f>, tt + 27rk 1 (s m + Ort^'y. 



+nkMs m + n m %)G mn (s n + n n %)\ ^ T ( 7 m ) (32) 

The only functions that are doubly (semi-)periodic are combinations of the theta functions 
(PP. After some algebra, we find that the set of functions 

/ \ / a+k 1 p+k 2 r \ 

%,r ( p J = e^ m ^0 f fata J (fciTlAnfeafi) (33) 

where p = 1,2, ...,k 2 and r = 1,2,...,^ with a, /3 G [0,1] solve the above algebraic 
conditions (0). Their inner product is given by 

(*Wi*W 2 )= / e" 2 ^ mG -^^Hv]/ p2ir2 ( 7 )|G|- i nd7 m d7 m (34) 

\G\ 2 3pi,p2^ri,r2 

The integrand is completely invariant under the translation Q28|), thus we restrict the inte- 
gration to one of the plaquettes P ( r y m = u m + Q ml vi with u, v G [0, 1]), the phase space of 
the 7's. 

Under a large gauge transformation 

Cst%r ( « ) ( 7 ) = e ^^^+^(^-»-/^)^ r (J V T) 

= E^W*^ ( ^ ) (7) (35) 

The matrix [C Sji ] PiP / forms a (/C2) 9 dimensional representation of the algebra ([3(]) of large 
gauge transformations. 

The parameters a and j3 appear as free parameters, but in fact they may be fixed such 
that we obtain a modular invariant wave function. The modular transformation (|5|) on our 
set of functions ([53|) is 

where 7' = (Cfi + £>)~ lT 7, fi' = (AH + £)(Cfi + D)" 1 and is a phase that will not 
concern us here (and G' lm = [(Cfl + -D) _1 ]z r G> s [(Cf2 + _D) -1 ] sm ). Most important are the 
new variables 

a' = Da -CP- ^(CD T ) d (3' = -Ba + A(3 - ^(AB T ) d 
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where (M) d mean [M] rfrf , the diagonal elements. 

A set of modular invariant wave functions JIJ |5j, U can exist only when k\k 2 is even, 
where we set a = [3 = (and also = 0). In the case of odd kik 2l we can set a, [3 to either 
or |, which amount to the addition of a spin structure on the wave functions. This will 
increase the number of functions by 4 9 which will now transform non trivially under modular 
transformations. 

Considering a set of point charges leads to the set of wave functions 



%,r f p J (f,7,W =exp 



rt _ 2ixi 
7rkj m lm + 2vr 7 m J o Q m - j m )dt' - — E 



+Yk ? - %(o)) + ^ I (j m - J m )df ■ jT (J™ - r)dt' 

( a+fcip+fc 2 r \ *f 
fc j| 2 j (A:i7 m -A; 2 ^ j^dt'l^feO) 



(36) 



The wave function depends on charge positions through the integrals over the topological 
components of the current j m ,j m , and through the function 9ij(t) — 6^(0). Consider for a 
moment motions of the particles which are closed curves, and are homologically trivial. 
We focus first on the integrals over j m ,j m . If, for example a single particle moves in a 
circle, we find that the integral of these topological currents vanishes, we conclude that these 
currents contribute nothing additional to the phase of the wave function under these kinds 
of motions. The function %(t) — %(0) must be treated differently here, because it has 
singularities when particles coincide, and thus, while motions that encircle no other particles 
may be easily integrated to get zero, this is not true when other particles are enclosed by 
one of the particle paths, and the result is nonzero in this case, in fact it is 2ir. Nevertheless, 
this function is still independent of the particular shape of the particle path. Actually the 
definition of in term of the prime form E(z,w) is just the generalization to an arbitrary 
Riemann surface of the well known angle function on the plane, that is as the angle of the 
line joining the particle % and j compare to a fixed axis of reference, determined by zq here. 
Thus, we may conclude that, under permutations of particles of charge q, the wave functions 
defined here acquire the phase a = e~£ q . 

For homologically nontrivial motions of a single particle on M. , the current integral 
Iaj l (t')dt' will in general change as ^j\t')dt' — > ^j l (t')dt' + s l + Q lm t m , where s l and 
t m are integer-valued vectors whose entries denote the number of windings of the particle 
around each homological cycle. However, now, for multi-particle non-braiding paths, there 
is no contribution comming from 9^. Thus, the wave functions become 



# p ,r(£) = exp 



27T7 27T? ? 

-rr m s m - —{{a - k 2 a ) m s m -{(3- k 2 (3 ) m t m ) - —J 
k k\ Ik 



■ * PlP+ t(0) = J2l B s,t}r,r>* P A°) ( 37 ) 
r> 

with J2i q% Izo (0) u l = a l + n lm f3 0m and where J = J2i (/*(*) - /<(0)) is a self-linking term. 
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The self-linking contribution in plays a very important role. To construct a framing 
of the curve z iy consider transporting a vector along the curve using a metric compatible 
connection on A4. If this operation were performed on an open manifold, on a curve C 
enclosing a region U, then the vector, after transport around the curve C, would suffer a 
deflection by an angle given by a line integral of the connection around this curve, which 
equals J v R, where R is the curvature of the connection on M.. However, on a compact 
Riemann surface this operation may equivalently be viewed as a transport backwards along 
the curve C, encircling the complement M. — U, with the result — j M _u R- These results are 
not the same, but differ by J M R = 47r(g — 1) = 2nx(g)- Our problem with this calculation 
is that the Riemann surface can only be properly defined using a set of patches, and viewed 
from different patches, parallel transport of a vector has different meanings. This ambiguity 
complicates further the framing ambiguity already mentioned for which, we recall, a special 
choice, a constant framing implying J = 0, made our theory into a topological field theory. 
However, this choice can be defined only on a single patch - thus, for a consistent topological 
field theory, we must ensure that the additional phase ambiguity which arises in the wave 
function is rendered harmless, that is, 

e ^G>-i) 9 ? = ! ( 38) 

These are fundamental constraints, for g ^ 1, on the set of possible charges depending on 
ki and g. 



The matrices fl37|) satisfy the cocycle relation 

B Sutl B S2 , 2 = e-^^-^B S2 , t2 B Sl>tl (39) 

(to be contrasted with the large gauge transformations cocycle (0)). Together with the 
permutation phase a, these matrices are the result of the action of elements of the (permuted) 
braid group on the particles which form the external sources in our theory, and we might 
reasonably expect them to be a representation of the braid group. In fact, let the integer 
vectors s l , t m denote vectors that are in all entries except for the Zth and mth, respectively, 
and 1 at the remaining position. Then with the identifications a\ = B s i , f3 m = B £ m , it 
is easy to check that we recover all of the conditions that appear in the work of Imbo and 
March-Russell ((H), that is their equation (2.6). We do not quite find their constraint (2.5), 
which restricts the value of the phase a as <r 2(n+9_1 ' ) = 1. From (|3"8D, the above constraint is 
reduced to cr 2n = 1. Since they are considering a set of identical particles, we would have to 
work in the sector of total charge Q = nq. Actually, this constraint is a global condition on 
the representation of the braid group. Here we have to suppose that there is an additional (or 
a group of) particle of charge — Q present. The work in this non-zero charge sector is under 
progress now and will be published in a later paper. We expect to recover this additional 
constraint. 

We have quantized Abelian Chern-Simons Theory coupled to arbitrary external sources 
on an arbitrary Riemann surface, and solved the theory. We find that the presence of 
nontrivial spatial topology introduces extra dimensionality to the Hilbert space separately 
for the large gauge transformations and the braid group. We find the constraints (|38|), 
relating the charges, k, and g such that we recover a consistent topological field theory 
representing the braid group on A4. 
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